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We study the Bogoliubov spectrum of an elongated Fermi superfluid confined in a one-dimensional 
superfluid along the Bose-Einstein-condensate (BEC)-Bardeen-Cooper-Schrieffer (BCS) crossover. 
We derive analytic expressions for the effective mass and the Bogoliubov excitation spectrum of 
the axial quasiparticles along the crossover based on the hydrodynamic theory. Our investigation 
reveal interesting signatures of BEC-BCS crossover in an optical lattice which deserve experimental 
investigation. 



I. INTRODUCTION 

The experimental realization [T, '2*1 of optical lattices for fcrmionic isotopes such as ^Li or '^^K is stimulating new 
perspectives in the study of superfluidity in these systems. An increase in the superfluid transition temperature when 
using potentials created by standing light waves has been predicted Q. It has been shown that for values of the 
Fermi energy above the first Bloch band, the center of mass motion of a Fermi gas trapped in an one-dimensional 
periodic potential is overdamped in the coUisional regime due to Umklapp process Q. A hydrodynamic theory of 
Fermi superfluids in the presence of optical lattice has been developed in the weakly interacting BCS limit 5] . 

In such a fermionic system, it should be possible to adjust the interaction strength and light intensity to tune the 
system continuously between two limits: a Bardeen-Cooper-Schrieffer (BCS) type superfluid (involving correlated 
atom pairs in momentum space) and a Bosc-Einstein condensate (BEC) in which spatially local pairs of atoms are 
bound together. This crossover between BCS-type superfluidity and the BEC limit for a dilute gas of fermionic atoms 
has been of recent theoretical interest In particular lot of theoretical attention has focussed on the collective 
excitations at the BEC-BCS crossover y\. The first experimental results on the collective frequencies of the lowest 
axial and radial breathing modes of ultracold gases of ^Li across the Feshbach resonance have also become available 
Q. In atomic Fermi gases, tunable strong interactions are produced using the Feshbach resonance Feshbach 
resonance, occurs when the energy of a quasibound molecular state becomes equal to the energy of two free atoms. 
The magnetic field dependence of the resonance allows precise tuning of the atom-atom interaction strength. Across 
the resonance the s-wave scattering length goes from large positive to large negative values. The fermionic system 
becomes molecular BEC for strongly repulsive interaction and transforms into a BCS superfluid when the interaction 
is attractive. Recent experiments have entered the crossover regime and yielded results of the interaction strength by 
the cloud size and expansion (lo| . 

In atomic Fermi gas experiments, BEC-BCS crossover regime in optical lattices has not yet been demonstrated. 
Considering the fact that experiments in this area are making rapid progress, we were motivated to study for the 
first time the low energy Bogoliubov spectrum of a cigar shaped superfluid Fermi gas confined in an one-dimensional 
optical lattice along the BEC-BCS crossover regime using the hydrodynamic approach. For fcrmions confined in a 
trapping potential, the density profile changes slowly in space if the particle number of the system is large enough. 
Under such conditions, a local density approximation can be applied to the state and a hydrodynamic approach can 
be adopted to investigate the low energy collective modes. Expressions for the effective mass and various collective 
modes are new results of this work. 



II. THE EFFECTIVE ACTION AND THE HYDRODYNAMIC EQUATIONS IN AN OPTICAL 

LATTICE 

We consider a cigar shaped dilute ultracold fermionic gas trapped in an one-dimensional optical lattice. The optical 
lattice is formed by two counterpropagating laser beams, for example in the z direction. 
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Here, d is the lattice period and s is the dimensionless amphtude of the lattice potential. Ej^ ~ 



2m(P 



is the recoil 



energy [ur — is the corresponding recoil frequency) of the lattice. In addition to the optical lattice, we also have 

TTl 

a harmonic trap V/,,o(?', ^) ~ "2" (^^^^ ^ '^z^^)- Mass of the fermionic atom is m. and lo^ are the radial and the 
axial trap frequencies. We take lOj. > ^z- The harmonic oscillator frequency corresponding to small motion about 

the minima of the optical lattice is Ug ~ — — . t^s >> so that the optical lattice dominates the harmonic 

potential along the z-direction and hence the harmonic potential is neglected. The strong laser intensity will give 
rise to an array of several quasi-two dimensional pancake shaped condensates. In writing down the effective action, 
we will follow Wouters et al. [llj. When an optical lattice is present along the z-direction, we can decouple the free 
motion in the a;,y-plane from the tunneling motion in the z-direction. The partition fimction for a system consisting 
of layers of 2D fermions is 



where the action is given by 



(2) 
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Here the three vector notation is used, r — {x,y,T). (3 = l/fc^T, T is the temperature. Also x = x^y. The 
field Tpj^a-if) belongs to a fermion of mass m in layer j and spin or |). U is the attractive strength between 



(I) j;i2]. Kxt(j) is the 



the fermions. The interlayer Josephson tunneling energy is J = sEb. — ij exp ^— a 

external potential acting on each layer and can be parabolic in addition to the optical lattice. In order to get rid 
of the quartic term in equation (3) one needs to perform a Hubbard-Stratonovic (HS) transformation and introduce 
the HS fields Aj(r) cx< V'jV'j >j after which the integration over fermionic variables is performed. Our goal 

is an investigation of the superfluid properties of the ultracold Fermi system using the hydrodynamic approach. The 

hydrodynamic interpretation of |A,(r)P is that it represents the density of fermion pairs, whereas v{r) = ^ j( ^ jg 

m 

interpreted as the superfluid velocity field. The Hubbard-Stratonovich transformation transforms equation (2) into: 



Z 



(4) 



with 



^ / dT I dx 



|A, (r) 



U 



a=±l \ m / 

-A, (r) ^pI^ (r) ^pI^ (r) - A] (r) t/;,-^ (r) 1 (^) 
In order to have information about the physical density, we multiply the partition function by constant C 



(5) 



C = J VQ {r)Vp, (r) exp | - ^ ^ dr 1 dx (r) 

X (r) - {r) V-j.T {r) - V-jj (r) V-j^i {r) | . 



(6) 
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Carrying out the functional integral over the new HS field Q (r) gives 5[pj (r) — ip] ■] {r) tpj^^ (r) — V"! x (^) i^jd (^)] 
and thus pj (r) corresponds to the physical density of the system along any path. In order to investigate the BCS 
gap and the phase we separate the complex field Aj(r) in a modulus and a phase Aj (r) = |Aj (r)| e*^^''"-' and also 
transform the fcrmion fields as V'i.o- (r) {f) e'^^'-''-'/^. Additionally, we shift the field iQ (r) according to 



iQj {r) — > iC,j (r) + ihBr- 



8m 



(7) 



and use the Nambu spinor notation rjj (r) = (^j,] {f) , V'j | (^) 

After this procedure and a further path integral over the fermion fields results in the partition function to be written 



j 2?|A,(r)|P0,(r)PO(r)Pp,(r)exp{-5effA}, 



(8) 



where 
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and 
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Here Ui are the Pauli matrices. Now to determine the low energy dynamics of the density and the phase of the 
superfluid, the paths along which 9j{r) and Pj{r) vary slowly in comparison to the fermionic frequencies (Fermi energy 
and binding energy) are of importance. Along these paths we have to make saddle point approximations for the fields 
|Aj(r')| and C,j{r) as: 



|A,(r) 



Af (.) 
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where 



Af ) (.) 



and C-*^^ (r) also vary slowly in comparison to the fcrmion frequencies. Since we are interested in the 



low-energy regime of Seff, corrections that go beyond this regime are neglected 115|. Since we are dealing with T = 
dynamics, the normal component is absent and the superfluid is stiff. This means that we neglect terms proportional 
V^6'(r) 

to , which is equivalent to V.t/ = (we consider steady flow). 

m , , , , 

Expanding the effective action (equation 9) and the Green's function (equation 11) around the saddle point values 



,(0) 



and g 



(0) 



leads to the saddle point effective action as a sum of contributions independent of J and 



tunnehng contributions(S'*y ) : 
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where, 
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The saddle point equations are: 

1 

U 



(fk l-2nF[Ej{k)] 
2Ej (k) ■ 



2nF [Ejik)] - 1 



{2Try \2m ' ' J {, E, (fc) 

with npiE) — l/(e^^ + 1) the Fermi-Dirac distribution function and Ej (fc) the local BCS energy defined by 



iQ '{r}] + A '{r) . 



2m 



(14) 



(15) 



(16) 
(17) 



(18) 



In the limit of T = and Thomas Fermi (TF) approximation (neglecting the kinetic energy), the local BCS energy 
is rewritten as 



E,^,lz- + A^. 



(19) 



Here, we have used the saddle point value in the absence of vortex |A°(r)| — A. The first saddle point equation 
(equation 16) corresponds to the BCS gap equation, whereas the second saddle point equation leads to the BCS 
equation fixing the Zj{r) in each layer at T = and TF approximation as: 



PjU A2 _ pjU 



(20) 



where E^^ is the binding energy of the molecule in each layer and we have defined the gap A as: 



A = 



E^PjU 



(21) 



This definition of the gap is consistent with the definition given in Ref. [13[ for small number of atoms per lattice 

A2 



site, if we identify the binding energy as E^^ ~ U. We also now define, fJ-eff = zj — fi = ^ ^ /i as the 



2Ei^ 



effective chemical potential in each layer. We now need to calculate the terms in the effective action that couple the 
different layers since we want to study the current perpendicular to the layers in which the atoms are confined. The 
tunneling contributions in the effective action S'^^'Y]i>j+i '^^'^ be treated perturbatively. The lowest order perturbative 
expansion of the action with respect to the tunneling part yields: 



Qtunnel 



I dr dx J j^j+ipj cos { 



(22) 
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where 



2EID ^ 

The binding energy of the molecule in each layer is given by [l 



El"" = 0.583 Er exp ( ] ■ (24) 



A is the wavelength of the laser light and is the scattering length. Having obtained the saddle point effective 
action S^j:j: that depends on 0j{r) and Pj{r), we can now derive the equations of motion for the variables dj{r) and 
Pj(r) from saddle point effective action through the extremum conditions dS^^j/dOj — and the number equation. 
This leads to the following hydrodynamic equations for Pj{r) and Oj{r): 



2~m^ (P^^^^) + -^J'J-i^J ^^"^(^J " ^J-i) " •^J+ij/'j sin(6'j+i - 0j), (25) 

h do fi^ 

Note that in equations (25) and (26), we have not explictly written the spatial dependence of pj and Oj since it is 
understood that they depend on r. In the next section our starting point will be the above hydrodynamic equations 
and solve the corresponding equations of motion for the density and velocity fluctuations. 

III. MULTIBRANCH BOGOLIUBOV SPECTRUM 

The equation of state enters through the density dependent chemical potential. We assume the power-law form 
of the equation of state as /ie//(/3) — Cp'^ Q. 7 is an effective polytropic index. The polytropic approximation has 
the advantage of allowing one to get analytical expressions for the eigenfunctions and eigenfrequencies of collective 
modes for various superfluid regimes in a unified way. At equilibrium, the two dimensional density profile takes the 

form at each layer pi^W) = i ^"/,^ \ (l — f^)^^'', where r = and R ~ 1/— ^— Linearizing around equilibrium. 

Pi — Po + Spj, 9j = 59 j and peffipj) — Peffipo) H — \p=p„Spj. The equations of motion for the density and phase 
fluctuations are 

' E^ 



E2D -r -b 

In deriving the above equations,we have assumed that the optical lattice is deep so that Jj±ij << Peff- The 
second order equation of motion for the density fluctuation is given by 



where 
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/2 

F = ^ (30) 

(#+Er) 



A further simplification yields 



We assume a normal mode solution of the form 

5pj{f,z,t) = Spo{f) exp i {u}a{k)t — jkd) . (32) 

This yields 

- ^I5p.{f) - ((1 - f^)^/W.(l - r-Y-''^'5p,{r)) l^^^lll sin^ (^^) (33) 

where, uj = — . Here a is a set of two quantum numbers: radial quantum number and the angular quantum 
number m. 

For fc = 0, it reduces to a two-dimensional eigenvalue problem and the solutions of it can be obtained analytically 
by generalizing the method of (l6j to any polytropic equation of state. The energy spectrum is given by 

= \m\ + 2nr hirir + \m\) + 1) (34) 
The corresponding normalized cigcnfunction is given by 

^ A (1 -f 2) (1/7-1) ~|™| pa/7-1, (2^:2 _ exp iim(j)) (35) 

where P°'''(x) is a Jacobi polynomial of order n and is the polar angle. Also, the normalization constant A is 
given by 



,2 22-2/7 [r(n,, + l)]2r(l/7)r(2/7 + 2n, + |m|) 



V^i?2 r(l/7-l/2)[r(l/7 + n0]2r(2n, + |m| + l)" ^ ' 

For fc ^ 0, we expand the density fluctuation as 

SpQ ='^baSpa{f,(j>). (37) 

a 

Substituting the above expansion into equation (33), we obtain 

= p2 _ ^ 2nr{-i{nr + \m\) + 1)] 

- Bosin2 (^^^]6„ + Bosin2 (^^^ ^M„„-6„'- (38) 



where 



= (39) 
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X P^7"'"''^(2f2_l)p^70,|m|)(2f2_i)^ (40) 

where 70 = I/7 — 1, = ttR^A^. Equation (40) is the central resuh of this work, which is derived for the first 
time. 

Before proceeding ahead, let us have a look at typical experimental values. Let us now consider the experiment 
of Greiner et al. [l7|- They had prepared ultracold gas of fermionic "^^K atoms in the lowest energy spin states 
1/ = 9/2, m/ = —7/2) and |/ = 9/2, m/ = —9/2), where / is the total atomic angular momentum and m/ the 
magnetic quantum number. The Feshbach resonance occurs at a magnetic field Bq = 202.1 ± 0.1 G and has a width 
of ui = 7.8 ± 0.6 G. When B < Bq, the s-wave scattering length as is positive and a weakly bound molecular BEG 
is formed. At B ~ Bq (unitary limit), ±00, which corresponds to a very small binding energy ( see equation 

24). Beyond the resonance {B > Bq), is negative and the binding energies in the BGS(-E^^),umtary(iJ,^^-) 
and the BEG(i?^p^) region are related as E^^^ < E^^^ < El^^. The scattering length Os is determined as = 

174ao ( 1 + 7 r ) 17], where ao is the Bohr radius. 

V [Bq^ - B)J ^ ^ 

Since we are interested in the crossover, we will focus our attention near the unitary limit. On the BEG side of the 
unitary limit, if we take B — 201.2 G < Bq, we find that the scattering length on the BEG side is a''f w 3.78 x 10~^m 
and E^^^ « 4.26Er. In the BEG region, the chemical potential is ^ = -E^^^/2, A » Ep and Ep < Ej^^. Using 

these values, we find that El^^ » On the BGS side near the unitary hmit, B = 202.2 G, a!f '' « -3.60 x 10^'^m 

^bec 

and El^^ « 0.22£;i^, = ^^^(Fermi energy), A^O,Ef » El^^. 

To calculate the multibranch Bogoliubov spectrum, we need to know the adiabatic index 7. 7 = 2/3 denotes the 
unitary as well as the BGS limit For the BEG side of the unitary limit 7 > 2/3 and for the BGS side of the unitary 
limit we take 0.6 < 7 < 2/3. 

In Figure [U we plot the lowest (n^ = 0, m = 0) mode for the BEG side of unitary limit (thick line, = 

Er 

1,7 = 0.8, A = 2EF,Elg = 8Ef) and BGS side of unitary limit (thin line,^:^ = 25E^^^,— = 1,7 = 0.6, A = 

Er 

0.4i5i?).Glearly, we find that in the entire Brillioun zone, the frequency of the phonon mode in the BGS side of unitary 
limit is greater than that in the BEG side of unitary limit. In the limit of long wavelength, the rir = mode is 

phonon like. On the BEG side ujbec ~ k^j — ir^—^, where rn^^^ = — '"^^ is the effective mass in the BEG side of 
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2.2 




3 

kd 



Figure 2: Plots of the monopole modes (m — 0,nr — 1) in the BEC regime (thick line,-^^ = 1,7 = 0.8, A — 2Ef, -Eti^ = iEp 

Er 



) and BCS regime (thin line, Ep = '^^El 
wavelength and the short wavelength behaviour. 



E^ 



1,7 = 0.6, A — OAEf)- One can clearly see the difference in the long 



unitary limit. On the BCS side the low wavelength behaviour is Wbc 



/7(2-7)^e// 



"bcs 



where m 



bcs 



In Ref. it was found that the effective mass increases as the density of the particle increases. This observation 
is in accordance with our result, if we substitute / El j^^ from equation (21) in the expression for m^^^. Note that 
the expression for mj^^ coincides with that found in and m^^^ > m^^^. By gradually changing the magnetic field 
when we go from the BEC regime to the BCS regime, we observe an increase in the phonon mode which is expected 
for a coUisionless Fermi gas, where the elastic collision rate is strongly reduced by Pauli blocking. In the absence of 
coupling between the axial modes with the density along the radial direction (for the rt^ = mode), the effect of Pauli 
blocking is strong and ujbec < ^bcs hi the entire Brillioun zone, the coupling between the axial modes and the radial 
density enhances the elastic collision rate. 

In Figure[2l we show the monopole modes (m = 0, = l)in the BEC side of unitary limit (thick line) and the BCS 
side of unitary limit (thin line). The monopole modes show a very peculiar behaviour. In the long wavelength region, 
LUbec > ^bcs- As we movc away from the center of the Brillioun zone, there is a cross over and ujbec < i^bcs- If we look 
at equation (21), we find that in the long wavelength region, the term 2nr['y(nr + 1)] (this part indicates the coupling 
between the radial density and the axial frequencies and is finite for ^ O)dominates over the part determined by 

the optical lattice (1 — Ma^a) Bq — - — . Consequently, the elastic damping rate is high and cj^ec > ^bcs as expected 

from experiments on elongated Fermi gases because the term 2nr['y{nr + 1)] is greater in the BEC region as 
compared to that in the BCS region (since "/bee > Jbes)- On moving away from the center of the Brillioun zone, the 

term proportional to sin — ) starts dominating and as we go from the BEC to BCS regime, we probably enter 



the coUisionless phase and the elastic collision rate decreases and as a result cobee < ^bes- Note that Pauli blocking 
reduces the binding energy because of which E'^^ < E^^. For both the phonon modes and the monopole modes, 
the frequency {u>uni) in the unitary limit(7 = 2/3) will lie between ujbee and uJbes- Eventhough 7 = 2/3 in both the 
unitary and BCS limit, oJmii < oJbes because from equation (21) we notice that E^^^ > i?^^. Our analysis of the 
low energy Bogoliubov modes in the BCS-BEC crossover was obtained in the Thomas-Fermi approximation, which 
is valid only in the large density (weak coupling) regime. The important excitations in the BCS limit involve broken 
pairs. With increaing attraction these are pushed to very high energies, and in the Bose limit it is the collective modes 
which are the dominant low-energy excitations. Thus a proper description of the intermediate coupling regime must 
include both broken pairs and the collective modes. It is clear that present approach is not suitable to describe the 
intermediate as well as well as the strong coupling regime. A better method is the LDA (local density approximation) 
which is known to give satisfactory results in the intermediate as well as the strong coupling regime jlSj . The spectra 
of the phonon and the monopole modes in the different regimes can be observed in the Bragg scattering experiments 
as these spectra have been observed in Ref. [l^ for weakly interacting BEC. By measuring the sound velocity in 
pulse propagation experiments and by observing the low-energy Bogoliubov spectrum in the Bragg spectroscopy, one 
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can make a clear identification of various superfluid regimes along the BCS-BEC crossover. The results presented in 
this work may be useful for guiding experiments, which look for signatures of BEC-BCS crossover in optical lattices. 



IV. CONCLUSIONS 

We have studied the Bogoliubov spectrum of an elongated Fermi superfluid confined in an one-dimensional superfluid 
along the Bose-Einstein-condensate (BEC)-Bardeen-Cooper-SchriefFer (BCS) crossover. Using the hydrodynamic 
approach, we have analytically calculated the effective mass and the multibranch Bogoliubov spectrum in the BEC 
and BCS side of the unitary limit. We have shown that the effective mass increases as the system crosses from the BCS 
side to the BEC side. The Bogoliubov axial excitation frequencies on either side of the unitary limit show a strong 
dependence on the coupling with the radial density and the binding energy and thus provide valuable information on 
the physical behavior of the system. The frequency of the phonon mode in the BCS side is greater than that in the 
BEC side. On the other hand near the center of the Brillioun zone, we show that the monopole frequency on the BEC 
side is greater than that on the BCS side but as we go towards the edge of the Brillioun zone, the monopole frequency 
on the BCS side becomes greater than that on the BEC side. The various Bogoliubov frequencies calculated here can 
be measured by Bragg scattering experiments. 
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